Entangled state engineering of vibrational modes in a multi-membrane optomechanical system 
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We propose a method to generate entangled states of the vibrational modes of N membranes which are cou- 
pled to a cavity mode via the radiation pressure. Using sideband excitations, we show that arbitrary entangled 
states of vibrational modes of different membranes can be produced in principle by sequentially applying a 
series of classical pulses with desired frequencies, phases and durations. As examples, we show how to syn- 
thesize several typical entangled states, for example, Bell states, NOON states, GHZ states and W states. The 
environmental effect, information leakage, and experimental feasibility are briefly discussed. Our proposal can 
also be applied to other experimental setups of optomechanical systems, in which many mechanical resonators 
are coupled to a common sing-mode cavity field via the radiation pressure. 
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I. INTRODUCTION 

Entanglement is one of the most important resources for 
quantum information processing [1]. Entanglement of the in- 
ternal freedoms of microscopic systems has already been pre- 
pared experimentally [2], for example, the polarizations of 
photons yj], the electronic states of atoms Jj] and spin states 
of ions J5[]. The entangled states in macroscopic supercon- 
ducting quantum systems have been experimentally demon- 
strated (0, 0] . Recently, the quantum properties of the me- 
chanically vibrational modes are extensively studied from 
external degrees of freedom of microscopic particles (e.g., 
trapped ions 18[]) to macroscopic objects ||9|]. 

The entanglement generation and arbitrary quantum state 
preparation of vibrational modes in microscopic systems 
(e.g., trapped ions) have been studied in both experiments and 
theories (e.g., see Refs. II 1 Ql - 4 1 311 ) . Thus we would like to 
know whether the entanglement can be generated in the sys- 
tems of macroscopic mechanical resonators. The research 
on the coupling between superconducting quantum devices 
and macroscopic mechanical resonators shows that arbitrary 
phonon states can be produced in principle by using a pro- 
posed method lfl4ll when a macroscopic mechanical resonator 
is coupled to a superconducting qubits 01511 . also the gen- 
eration of the squeezed and entangled states of two vibra- 
tional modes has been proposed by coupling two macroscopic 
mechanical resonators to the superconducting quantum de- 
vices llql . However, their experimental realizations are still 
very challenge. The main obstacle is whether the macroscopic 
mechanical resonator can be in its ground state. 

The ground state cooling IU7I432I1 of the macroscopic me- 
chanical resonators has been theoretically studied and exper- 
imentally demonstrated by coupling them to a cavity field 
via the radiation pressure. Thus the optomichanical systems 
(see, reviews j33ll34lo provide a very good platform to study 
quantum mechanics at the macroscopic scale. The stationary 
entanglement between mechanical and cavity modes in the 



* yuxiliu@tsinghua.edu.cn 



optomechanical systems has been studied JUrHH], and such 
continuous variable entan glem ent can be used to implement 
quantum teleportation 1 40-42]. Also both tripartite and bipar- 
tite entanglement between mechanical modes and other de- 
grees of freedom can be generated in the optomechanical sys- 
tems I43l449ll or the hybrid optomechanical systems with an 
atomic ensemble II50I - I59I1 or a single-atom Il60l 16111 inside the 
cavity. Moreover, the cavity field mediated entanglement be- 
tween two macroscopic mechanical resonators in the steady 
state ll62U72ll has also been theoretically explored. How- 
ever, the coherent engineering of arbitrarily entangled phonon 
states of the macroscopic mechanical resonators is still an 
open question. 

We have studied a deterministic method, which is differ- 
ent from the proposals on measurement-based on phonon 
state generation 117 3l to synthesize arbitrary nonclassical 
phonon states in optomechanical systems lf75ll . Recent stud- 
ies show that many macroscopic mechanical resonators can 
be coupled to a common single-mode cavity field ll76l478l[80l - 
H3l via the radiation pressure, which has been theoretically 
studied for selected entanglement generation llsill . synchro- 
nization |82ll and mechanical analogue of nonlinear quantum 
optics [83] of many mechanical modes, and also experimen- 
tally demonstrated tripartite mixing for one cavity mode and 
two mechanical mode in the system that a microwave cavity is 
coupled to two and more mechanical resonators fl84h . In such 
system, the single-mode cavity field is a very good candidate 
to act as a data bus for information transfer from one mechan- 
ical mode to another one. Motivated by those researches, we 
use a multiple membrane optomechanical system J76l - f78[ fsotl 
as an example to study entangled phonon state engineering, in 
optomechanical systems with the many mechanical resonators 
coupled to a single-mode cavity field, by using the sideband 
excitations and and the single-photon effect 1 79] induced by 
the photon blockade. In particular, we will study detailed 
steps on engineering several typical entangled phonon states. 

The paper is organized as follows. In Sec. II, the theoreti- 
cal model of the optomechanical system with multiple mem- 
branes inside a single-mode cavity is introduced. In Sec. Ill, 
we study a method to generate entangled states for the system 
parameters within so-called Lamb-Dicke approximation. As 
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examples, we show how to generate Bell, NOON, GHZ and 
W states. In Sec. IV, we study entangled state engineering 
beyond the Lamb-Dicke approximation for the strong single- 
photon optomechanical coupling. Finally, discussions on ex- 
perimental feasibility and conclusions are given in Sec. V. 



n. THEORETICAL MODEL 
A. Mode equations and transfer matrix theory 

The mode equations of optomechanical systems with one 
and two membranes inside cavity have been studied by using 
the boundary conditions and the Helmholtz equations [64, 76, 
l77tl . However, the mode equations become harder and harder 
to be solved with the increase of the membrane number. The 
transfer-matrix method has been used widely in optics to an- 
alyze the propagation of electromagnetic fields especially in 
multi-layer structures [85]. It has been used to study the scat- 
tering problems in optomechanical systems [86-88]. For the 
completeness of the paper, we will first derive intrinsic mode 
equations of the optomechanical systems with N membranes 
inside a cavity using the transfer-matrix method. 

As a schematic diagram in Fig.Q] we focus on an optome- 
chanical system with a cavity containing N non-absorptive 
membranes, that each has reflection coefficient R, mass Mi, 
position qi, and vibration frequency (i = 1, • • • , iV). We 
assume that the thickness of each vibrational membrane is 
much smaller than the wave length of the cavity mode, so 
the electric susceptibility in the cavity can be approximatively 
described by a sum of the dielectric permittivity with 5 func- 
tions iMB 



e (x) = e 1 + 



c N 

i=l 



5(x- qi 



(1) 



where e is the dielectric permittivity of vacuum, £ = 
I^J Rj (1 — R), and k = ui/cis the wave vector of the elec- 
tric field with the mode frequency ui and the speed c of light 
in vacuum. 

It is well known that the transfer matrix describing the elec- 
tric field through the empty space of length I is given by ll85ll 



M {k, I) 



cos kl r sin kl 

k 

— k sin kl cos kl 



(2) 



Let us now study the transfer matrix of the whole system by 
exploring the properties of the electric fields cross a mem- 
brane. The continuity of E (x) at the position of the ilh mem- 
brane (e.g. x — qi) is give as 



E 



(3) 



where E [q[) and E [qt) are the notations of the left- and 
right-hand limits of E (x) when x approaches qi. Using 
the Helmholtz equations, E" (x) — — w 2 ^o£ (x) E (x), the 
derivative relations of E(x) in the left- and right-hand of the 
ith membrane at the position qi is given as 



E' (qf)=E' (qr)-kCE( qi ). 



(4) 
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FIG. 1 : Schematic diagram of an optomechanical system with TV me- 
chanical membranes inside the cavity which is driven by a classical 
field. Here, the two mirrors of the cavity are fixed, go and gjv+i de- 
note the positions of the cavity mirrors, gi, • • • , (ft, • • • ,qN denote 
the positions of the N membranes. 



with E' (x) = dE (x) jdx. Eqs. (0 and © can be written in 
a matrix form as 



where 



^o=(_i cl ) 



(6) 



is the transfer matrix of the electric field through the zth mem- 
brane. Using Eq. (f2]i and (O, the relation of the electric fields 
at the left and right mirrors can be given by 



E(q N+1 ) \_ Y f E(q ) 
E'(q N+ i) - XN {E'(q ) 



(7) 



with the transfer matrix 



Xjy — 



Xn X12 
X21 X 2 2 



(8) 



= [] [M (k, q t+1 - q t ) Q (k, 0} M (fc, qi - q ) . 

i=l 

Eq. (O can be rewritten as 

E(q N+1 ) = x n E (q Q ) + X 12 E' (q ) , (9) 
E' (q N+1 ) = x 2 iE ( 9o ) + X22E' (go) ■ (10) 

If we assume that the electric field satisfies the standing wave 
boundary conditions 



E(q o ) = E(q N+1 ) = 0, 
then we obtain the intrinsic mode equation 

X12 (k,qo,qi, ■ ■ ■ ,qN+i) = 0. 
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FIG. 2: The frequency shift of the intrinsic optical modes (uj — 
w^)/wf for two membranes in the cavity as function of displace- 
ments of the membranes, (51 — ?i)/A n and (52 — (?2)/An, where 
A n = 2L/n, uji = ttc/L, U)° = nnc/L, q- A ~ q° = q° — qi = 
qi — go = L. The parameters are R — 0.7, n — 10 5 . 



functions of the displacements of the membranes when there 
are two vibrating membranes inside the cavity. Analytically, 
the frequencies of the cavity modes can be given as 



N 



u(Ui}) = «({«?}) +£* (1) (*-«?) 



JV 



(2) I 

Tij [<H 



under the condition (qi — q®)/\ <C 1, where 



(i) 
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(2) 



doj({qt}) 



{«=«?} 



dqidqj 



(16) 



(17) 
(18) 



A is the wave length of optical mode, and g° (i = 1, 
is the position of the ith membrane when there is no radiation 
pressure. 



We further study the intrinsic mode equation in Eq. (fTzt by 
several concrete examples. For N = 1, the wave vector k 
obeys the intrinsic mode equation 



CVo + CCi,i = 0, 



(13) 



with 



Ci, = sin A: (g 2 - go) , 

C1.1 = sink(qi -q 2 )smk(qi - q ) . 

For N = 2, the intrinsic mode equation for k is 



C2,0 + CC2,1 + C C*2 



0. 



(14) 



with 



C2,0 
C2,l 

C2.2 



sinfc (q 3 - q Q ) , 

sin k (qi - q ) sin k (qi - q 3 ) , 

+ sinfc (g 2 - qo) smk(q 2 - q 3 ) , 

sinfc(yi ~q )smk(q 2 - yi) sin fc (q 3 - q 2 ) ■ 



If there are N membranes inside the cavity, the wave vector k 
obeys the following equation 



JV 



0, 



(15) 



i=0 



where Cjv.i are functions of the qi, i = 0, 1, 2, • • • , TV + 1. 

The frequency dependence of the cavity modes uj ({qi}) on 
the positions of the N membranes can be obtained by solving 
the intrinsic mode equation in Eq. (fT2] > numerically or analyt- 
ically by using perturbation theory. Here, {qi} is an abbrevi- 
ation of {qi, ■ ■ ■ , <Zjv} for tne set of the positions of 
the membranes. For example, in Fig. [2] several intrinsic mode 
frequencies have been numerically simulated and plotted as 



B. Hamiltonian of the system 

Based on above discussions, the Hamiltonian of the op- 
tomechanical system with N membranes inside the cavity can 
be written as B91I1 

N r 2 1 

H = tw({q i })aU + Y J ^+2^^-^'' 

? — 1 1 

(19) 

where a (a 1 *) is the annihilation (creation) operator of the cav- 
ity field, pi is the momentum of the zth vibrational membrane. 
Substituting Eq. ( fT6b into Eq. ( fT9b , we have 



N r 



Pi 

2M t 



AT 



(20) 



! = 1 

JV 



The third and fourth term are the linear and quadratic interac- 
tions between the cavity mode and the vibrational modes of 
the membranes. 

Hereafter, we only consider that the frequency shift of the 
cavity mode is linearly dependent on the membranes' dis- 
placements. We also assume that the cavity is driven by an 
external field with the frequency uud and the phase <fid- Thus, 
we have the Hamiltonian of the driven system as below 

TV N 

Hd = Tkjj a a) a + ^^hbjjb\bi + ha^a cji (b\ + b^j 



i=l 



i=l 



(21) 
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il denotes the Rabi frequency of the driven field. We assume 
that both the frequency u d and the phase <fi d are controllable 
parameters such that they can be chosen as different values in 
the steps of the state preparation described below. For the 
simplicity, the frequency of the cavity mode is denoted as 
uj a = oj ( { qf } ) , and the coupling strength between the cavity 
field and the ith membrane is simply written as The oper- 
ators of the membranes are rewritten by the annihilation and 
creation operators bi = -J MiUH / 2h(qi — q®)+ipi/ \ / 2hMiU>i 
and b\ = \J MiUii/2h(qi — g?) — ipij 1 \J2hMiiOi. We now 
applying a unitary transformation to Eq. (fJT) 



JV 



U = cxp I q^a^~~^ 



i=l 



9 - ( b\ b, 



then the Hamiltonian in Eq. ( f2TT > becomes into 



JV 



H e g — h(uj a — A a^ a)a) a + TjS^u>jb\bi 



Mlo) 



exp 



where Ao = Ylz=ii9i /^i) characterizes the nonlinearity of 
the cavity field induced by the vibrational membranes, thus 
the nonlinearity increases with the increase of the number of 
the membranes. We call rji — gi/uJi as the Lamb-Dicke pa- 
rameter in analog to the trapped ions [10]. Both the strong op- 
tomechanical coupling and many mechanical resonator make 
the nonlinear term Ao guarantee the photon blockade, in this 
case, the driving field can be assumed to be coupled to two 
lowest energy levels |0) and |1) of the cavity field. Thus, the 
cavity field is confined to the states |0) and 1 1), and the Hamil- 
tonian in Eq. d23l becomes into 



' JV 

y^J}j (b\ -fy) -i (uidt + 4> d ) 



(22) 



(23) 



H.c. 



#two = h^Y<T z + h^LUiblbi 



N 



Mla + exp 



i=l 
' N 



E Vi (b\ - bi) -i (oj d t + <t>d) 



(24) 



H.c. 



by using the operators a z = |1) (1| — |0) (0| and a + = |1) (0|. 
Here, ujq = u) a — Ao. Hereafter, we denote the photon number 
states as |1) = |e), |0) = \g). The Hamiltonian in Eq. d24"i i 
can be further written as 



with 
and 



Htwo — Hq + H u 



N 

H = h-^-Vz + E fwibik, 

i=l 



(25) 



(26) 



JV 



xb\H? +H.c. 



3iM=0 



(27) 



From Eq. d27l i. we find that — ji \ phonons can be created 
(hi > ji) or annihilated (fcj < ji) from the ith membrane 
when one photon is annihilated in the cavity with the assis- 
tance of the external field. Below, we will show entangled 
state engineering for different vibrational modes of the mem- 
branes for two cases with or without Lamb-Dicke approxima- 
tion. 



III. ENGINEERING ENTANGLED STATES WITH 
LAMB-DICKE APPROXIMATION 

We first study the entanglement engineering under the 
Lamb-Dicke approximation condition gi /a;, <C 1 as for the 
trapped ion case [10], thus the Hamiltonian in Eq. ( |27] i can be 
written as 



JV 

i+E^( 6 M*) 



H.c. 



(28) 

up to the first order of r\i . In the interaction picture with V = 

exp (iH t/h) Hi n t exp(—iH t/h), we have 



V = ma+e-Wd 
+H.c, 



N 



e" iA ^ - he 



-iA*t 



i=l 



(29) 



with A c = uj d - u Q , A| = uj d ~ w - Ui, and A* = u d - 
ljo + uji. If the system satisfies the resonant condition either 
A c = or A^ = or A* = 0, and also the driving field is not 
very strong, then we have 

{cr+e-'^+H.c. w d = w , 
T/iCT+ftte-^i + H.c. cu d = luq + Wi, (30) 
r)icr+bie~"fir + H.c. uj d = oj - w,. 

with the rotating wave approximation. For convenience, the 
minus sign before rjia+bie~ 1 ^^ is absorbed by the phase <j>\. . 



A. The time evolution operators 

From the Schrodinger equation, the wave function of the 
system at the time t can be written as 



\^(t)) = U(t)\i/ J (0)), 



(31) 



where U (t) = exp (—iVt/h) is the time evolution opera- 
tor 111311 . By using the completeness relation 



\-oo e 



£ 52\s,{m j }){s,{m j }\=I, (32) 

{ mj }=0 s=g 



the time evolution operator can be written as 

-f-oo e 



U ® = E E^W^'WXa.WI, (33) 

{m 3 }=0 s=g 
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where |{mj}) is an abbreviation of the state |mi) <8> • • • ® 
® • • • <8> |wjv) = | mi , • • • , to at) for iV membranes. Here 
\s, {to.,}) implies that the cavity field is in the state s (s = e or 
s = g) and there are mj phonons in the jth membranes, {nij } 
denotes a number series, that is, {rrij} = mi,m2, • • • , to at. 

If the frequency of the driving field is resonant with the 
red-sideband excitation corresponding to the frequency of the 
ith membrane, i.e., u>d — — then the Hamiltonian in 
Eq. (l30b becomes 

V% = hn^a+he- 1 ^ + H.c.. (34) 
In this case, the time evolution operator is given as 



U (IMXMIW (35) 

m,i=0 { mj }=0 



where 

UZ(t) = cos(fi^.t) m 4 ) (fl,m 2 | 



sin 



K<) i 



1) <ff,'< 



+ cos (fi^i+i*) l e ' m «) ( e ' m *l 

-ie**' sin {^l nz+1 t) \g,mi + 1) (e,mi|(36) 

with the Rabi frequencies 

f4. = n^y/WTt, 0^. +1 = Qr?iVmi + 1. (37) 

When the frequency of the driving field is resonant with 
the blue-sideband excitation corresponding to the frequency 
of the ith membrane, i.e., u>d = + Wj, the Hamiltonian in 
Eq. (f30b becomes 

V*£ = fflma+ble-^ + H.c. (38) 
The time evolution operator of the blue-sideband excitation is 

u%(t)=V&%(t) v awxwi)« ( 39 ) 



m«=0 



{ mi }=0 



with 



-ie~^* sin (O l mi+1 t) |e, m 4 + 1) (ff, m,| 

+ cos (ST mi i) |e, m*) (e, mj| (e, m,| 

-ie^» sin (Q l ro .i) | 5 , m, - 1) (e, m,| . (40) 

Finally, if the cavity is driven by the classical field with the 
frequency cod = u>o, then the carrier process is switched on, 
and the Hamiltonian in Eq. ( f30b is given as 

V c = Mla+e-^ + H.c. (41) 

The time evolution operator of the carrier process is given as 

U c (t) = U c (t) Y, IMXMI (42) 



and 

U c (t) = [cos (fit) |ff) - ie~ 1 ^ sin (Sit) \e)] (g\ 

+ [cos (fit) |e) - ie 1 ^ sin (fit) | 5 )] (e| . (43) 

Based on above three different processes, we can in princi- 
ple engineer any kind of entangled phonon states. However, 
below we will only show the engineering of several typical 
entangled phonon states by controlling the evolution time and 
the frequency of the classical driving field. 



B. Generation of Bell and NOON states of two mechanical 
modes 

In this section, we study the generation of entangled phonon 
states 



\<p) = -L(\N,0) + \0,N)), 



(44) 



of two vibrational modes when there are two membranes in- 
side the cavity. Here \<p) denotes the Bell state when N = 1. 
However \ip) with N > 2 represents the NOON state which 
plays an important role in quantum metrology 09211 . |iV, 0) 
means N phonons in mode one and zero phonon in mode 
two. Our state generation below starts from the initial state 
H 1 (to)} = Iff, 0, 0) of the whole system. 

We now show how to generate a Bell state. First, a driving 
field satisfying the carrier process is applied to the cavity, then 
after the interaction time Atj = 7r/20, the system state at the 
time t\ = t + At! is 



|^(ii)) = |e,0,0) 



(45) 



here a global phase has been neglected. Second, the frequency 
of the driving field is turned to the red sideband corresponding 
to the frequency of the first membrane, i.e., ojd = (jJq — 0Ji- 
With an evolution time A^2, the system evolves to 

IV> (*a)> = (l - |Ci, | 2 ) V2 |e, 0, 0) + d,o Iff, 1, 0) (46) 
at the time t<i = t\ + At2 with the parameter 
Ci >0 = -ie^ sin(n\At 2 ) . 



(47) 



If the time duration and the phase of the driving field are cho- 
sen as At2 — 7r/4fi}, phase (j>\ = n/2, then Ci.o = 1/V2, 
and the state of the system becomes into 



1^*2)} = ^|e,0,0) + i= |ff, 1,0). 



(48) 



In the third step, the frequency of the driving field is tuned to 
the red sideband corresponding to the frequency of the second 
membrane, i.e., ujd = — ^2- With the evolution time At%, 
the system evolves into 



m,)) = ± 



1 

71 



l-|Co,i|' 
9,1,0) 



1/2 



|e,0,0>+ Ceils, 0,1 

(49) 
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FIG. 3: Schematic diagrams for preparing NOON states 
(|2, 0) + |0, 2}) /y/2. The letters 71, B and C represent the process 
that the driving field is tuned to the red-sideband, blue-sideband or 
carrier processes, respectively. Each number of the subscript denotes 
the sideband excitations corresponding to the frequency of the ith 
membrane. 



at the time t% = t?. + At% with the parameter 



Co, 



-re 



(n?At 3 ) 



(50) 



(I) N is an even number 



|g,o,o; 



RilR^j)!*" 2 -" R, 
„|g,N/2.N/2> ► |c,N-l,0> >-|g,N,0> 



► |e,N/2-l,N/2> 



► |g,0,N> ►|g,0,N> 

R,(R,R,) IN/:M > R, 



(II) N is an odd number 



|g,o,o; 



-|e,(N-l)/2,(N-l)/2> ■ 



At, 



(R,Ri) (N_ly2 R, 
.|e,(N-l)/2,(N-l)/2> ► |e,N-l,0> ►|g,N,0> 
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FIG. 4: Schematic diagram for preparing NOON states when (I) N 
is an even number and (II) N is an odd number. The numbers in 
the superscript denote the times of the red sideband excitation (R), 
blue sideband excitation (B) and carrier process (C). The subscript 
denotes the sideband excitations corresponding to the frequency of 
the ith membrane. 



If we choose the time duration At 3 = tt/2Vi\ and phase = 
7r/2, then the state of the system becomes 



|^(t3)> = |ff>®-^(|0,l) + |l,0>). 



(51) 



Thus the system is deterministic ally prepared to a product 
state of the Bell state of two mechanical resonators and the 
ground state \g) of the cavity field. 

We now describe a detailed steps for preparing an arbi- 
trary NOON state of two vibrational membranes by the sim- 
plest example with N = 2. First, a driving field is applied 
to the cavity with the blue sideband excitation corresponding 
to frequency of the second membrane, i.e. ojd = wo + <^2, 
with an evolution time Ati = 7r/2ilf and choosing the phase 



= 37r/2, the system evolves to 

|V (ti)) = |e,0,l) 



(52) 



at the time t\ = to+Atx. Second, the frequency of the driving 
field is tuned to the red sideband excitation corresponding to 
the frequency of the first membrane, i.e., u)d =u)q — ljx. With 
an evolution time At2 = 7r/4fi \ and choosing the phase <f>l = 
7r/2, the state of the system becomes 



|V (t 2 )> = 4§(l2*M} + \e, 0,1)). 



(53) 



at the time t% = t\ + At%. Third, the frequency of the driving 
field is tuned to the red sideband excitation corresponding to 
the frequency of the second membrane, i.e., ujd = — <^2> 
with the time duration At$ satisfying 



sin(n?At 3 ) = ±l,sin(P^A£ 3 ) = ±1, 



(54) 



in the infinite approximation for = v2, which is an 

irrational number (see the Appendix B), the state of the whole 
system becomes 



hKts)) 



i 

71 



(|e, 1,0) + |s, 0,2)), 



(55) 



at the time t% = t2 + A^. In the last step, the driving field 
is tuned to the red sideband excitation corresponding to the 
frequency of the first membrane, i.e., ujd = — With an 
evolution time At4 = ir/2il}, and the phase <fi), = it/2, the 
state of the system becomes 



1 



|V(i4)) = |<?)®-/|(|2,0) + |0,2)), 



(56) 



at the time t 4 = t 3 + At 4 . Thus the NOON state for N — 2 is 
generated with the cavity field in its ground state \g). We have 
given a schematic diagram to summary all steps for generating 
the NOON state with N = 2 in Fig. [3] 

By using the similar steps as for generating the NOON state 
with N = 2, an arbitrary NOON state of two mechanical res- 
onators can also be generated by using the steps schematically 
shown in Fig. [4] The detailed steps are summarized as below. 
In the step (i), by alternatively applying a series of the red- 
sideband excitations and carrier processes, we can prepare the 
state 



IV(*l)) : 

for the even number TV or 

ltK*i)) = 



N 



N 



N - 1 7V-1 



e, ■ 



(57) 



(58) 



for the odd number N. 

In the step (ii), for the even number N, the system is driven 
with the red-sideband process R\ for the membrane one with 
the time duration At2 = 7r/4fiL, 2 , then we have 



1 



N N 
9l ~2'~2 



N N 

e 'T _1 'T 



(59) 



For the odd number N, the system is driven with the red side- 
band excitation i?2 for the membrane two with the time dura- 
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10 N 15 



FIG. 5: Fidelity for preparing NOON state (\N, 0) + |0, TV)) /s/2 
from TV = 2 to N = 20. 



\g, 0, 0, 0). The definitions of the GHZ state |V>) GHZ and W 
state \i/j) w are 



|V) GHZ = -/=(|o,o,o) + |i,i,i)). 



(63) 



and 



|V>) W = (|1, 0, 0) + |0, 1, 0) + |0, 0, 1)) . (64) 

Let us first show the steps for preparing the W state. In 
the step (i), the cavity field is driven by a laser field with the 
carrier frequency. With an evolution time Atj = 7r/2i7 and 
choosing the phase <fi c = 3tt/2, the state of the system be- 



|V(*i)) = |e, 0,0,0). 



(65) 



tion At 2 = 71-/402^ 



i)/2- 



then we have 



l^(*2)) 



1 

71 



N-l N - 1 



g, 



N - 1 N + l 



(60) 

In the step (iii), for the even number N, we alternatively 
apply the driving field for N/2 and (N/2) — 1 red sideband 
excitations for membrane two and one with the appropriate 
time durations, respectively, that is, an operation operator 
R^RiRz) 1 *' 2 ^ 1 is acted on the state \^p(t 2 )), then we have 



|^ 3 )) = -^(|e,iV-l,0> + |<?,0,iV>) 



(61) 



For the odd number AT, (N— l)/2red sideband excitations R 2 
and Ri are applied to the membrane two and one, respectively, 
then we also obtain the state shown in Eq. doTl ). 

In the step (iv), the system is driven with the red sideband 
excitation R\ for the membrane one with the time duration 
Ai4 = ir/2fl 1 N , we have 



\i/j(t i )) = \g)®-j=(\N,Q) + \Q,N)). 



(62) 



Thus the NOON state of the two mechanical resonators is pre- 
pared with the cavity field in its ground state \g). We note that 
there is information leakage in the third step because the times 
to the two states |e, N — 1,0) and \g, 0, N) are not synchro- 
nized. The fidelity of prepared NOON states from A^ = 2 to 
N = 20 due to such information leakage is given in Fig. [5] 
which shows that some states cannot be prepared with the 
hundred percent. The detailed discussion about this type of 
information leakage is given in the Appendix B. 



C. Generating GHZ and W states of three mechanical modes 



In the step (ii), the frequency of the driving field is tuned to 
the red sideband excitation corresponding to frequency of the 
membrane one such that ujd = wq— u\. With the time duration 
At 2 = [arcsin (l/V3)]/Oj and the phased = 7r/2,the state 
of the system evolves to 



IV>(* 2 )> 



1 



V3 



\g, 1,0,0) + W-|e, 0,0,0). (66) 



In the step (iii), the frequency of the driving field is tuned 
to the red sideband excitation corresponding to the mem- 
brane two such that uid = wo — w 2- With the time duration 
A^3 = 7r/4f^ and the phase $ = ir/2, the state of the sys- 
tem becomes 



IV>(*3)) 



1 



V3 



1,0,0) + |<7, 0,l,0) + |e, 0,0,0)). 



(67) 

In the step (iv), the frequency of the driving field is tuned to 
the red sideband excitation corresponding to the membrane 
three such that ujd = o->o — ^3- With the time duration 
At4 = 7r/20f and the phase 0j! = it/2, the state of the system 
evolves to 

|^(t 4 )) = | 5 )®^|(|l,0,0) + |0,l,0) + |0,0,l)). (68) 

Thus the system is prepared in a product state of the W state of 
three mechanical modes and the ground state \g) of the cavity 
field. 

Now, we are going to describe the steps of preparing the 
GHZ state, as schematically shown in Fig. [6] In the step (i), 
the cavity is driven with the carrier frequency. With the evolu- 
tion time Aii = 7r/20 and the phase <fi c = 3tt/2, the system 
evolves to 



\iHh)) = \e, 0,0,0). 



(69) 



We now turn to study the generation of the Greenberger- 
Horne-Zeilinger (GHZ) and W states i93ll of three vibrational 
membranes inside the cavity with the initial state \ip (to)} = 



In the step (ii), the frequency of the driving field is tuned 
to the red sideband excitation corresponding to the mem- 
brane one such that ujd = ujq — With the time duration 
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,1,0,0>— >- |g,0,0,0>— V |g,0,0,0>- 



c |->.|e,l,0,0>->- |g,0 
1,0,0> 

L>|g,I,0.0>->-|e,2. 



, At, At, 

t 3 t4 



|g,0,0.0>%. |g,0,0,0> 

.— >.|e,l,0,l>— ►fcl.l,l> 
Ri R : 

~At 6 ' At^ ' * 
t& t 7 Time 



FIG. 6: Schematic diagram for preparing GHZ state 
(|0, 0, 0) + |1, 1, 1)) /\/2 of three mechanical resonators. The 
numbers in the superscript denote the times of the red sideband 
excitation (R), blue sideband excitation (B) and carrier process (C). 
The subscript denotes the sideband excitations corresponding to the 
frequency of the ith membrane. 



At2 = tt/2H\ and the phase <\>\ — tt/2, the state of the sys- 
tem evolves to 



|V (ta)) = |5,1,0,0). 



(70) 



In the step (iii), the cavity field is driven again for the carrier 
process such that Ud = wo- Then, with the time duration 
A^3 = 71"/ 4£1 and the phase <\> c = 3tt/2, the state of the system 
evolves to 

|V (t 3 )) = -L \ gj 1, 0, 0) + i= |e, 1, 0, 0) . (71) 

In the step (iv), the frequency of the driving field is tuned 
to the blue sideband excitation corresponding to the frequency 
of the membrane one such that ujd = ujq + oj-y . For the time 
duration At.4 satisfying the condition 



sin(fj}A£ 4 ) = ±l,wa.(£l\M A ) = ±1, 



(72) 



with the infinite approximation (see the Appendix B), the state 
of the system evolves to 

|V (h)} = Iff, 0, 0, 0} + |e, 2, 0, 0) . (73) 

In the step (v), the driving field to tuned to the red side- 
band excitation corresponding to the membrane three such 
that ujd = ujq — ui3- With the time duration Atr, = 7r/2f2f 
and the phase <fil = tt/2, the state of the system evolves to 

|V (* B )> = A= IS, 0, 0, 0) + i= \g, 2, 0, 1) . (74) 

In the step (vi), the driving field to tuned the red sideband 
excitation corresponding to the membrane one such that ujd = 
too — uji. With the time duration /S.t§ = tt/2H\ and the phase 



37r/2, the state of the system evolves to 

\tl> (t 6 )) = -J= \g, 0,0,0} + -^=|e, 1,0,1) 



In the step (vii), the driving field is tuned to the red sideband 
excitation corresponding to the membrane two such that ujd = 
ujq — ll>2- With the time duration A<7 = 7r/2fq and the phase 
<f)\ = tt/2, the state of the system evolves to 



|V> (i 7 )> = 10)8-^(10,0,0) + |1,1,1)). (76) 



Thus the whole system is prepared to a product state of the 
GHZ state of three mechanical resonators and the ground state 
\g) of the cavity field. 

In principle, our method can be generalized to generation 
of the W and GHZ states of N membranes by sequentially 
applying a series of red-sideband excitations, blue-sideband 
excitations and carrier process with well chosen the time in- 
tervals, phases of the driving field. We have shown the steps 
for the generation of those states in the appendix A. 



IV. PREPARATION OF ENTANGLED STATES BEYOND 
THE LAMB-DICKE APPROXIMATION 

In the above, we have studied a detailed method of gener- 
ating entangled states in the Lamb-Dicke approximation with 
rji <C 1. However, in some optomechanical systems, e.g., a 
Bose-Einstein condensate served as the mechanical oscillator 
coupled to the cavity field Ir94l r95ll . the condition 77^ <§C 1 is 
not satisfied, and also with the experimental progress, the pa- 
rameter rji outside the Lamb-Dicke regime is possible to be 
achieved in the near future in the other types of optomechan- 
ical systems. For the parameter rji outside the Lamb-Dicke 
approximation, the higher-order powers of the Lamb-Dicke 
parameter should be taken into account. We now show how to 
generate entangled states beyond the Lamb-Dicke approxima- 
tion by using a similar but not same method given in Ref . Il2ll . 
Because the generation of the Bell state and the W state can 
use the same method as in the Lamb-Dicke approximation de- 
scribed in Sec. Ill, thus below we focus on the generation of 
the NOON state and the GHZ state. 

Beyond the Lamb-Dicke regime by assuming foj = ji + rii, 
the Hamiltonian in Eq. ( |27| ) can be rewritten as 

JV 

flint = hna+e-^+^Yle'^Hi + H.c, (77) 

i=i 



with 



+00 



+00 



Hi 



E E 

'»=— 00 j' j =max[0,' 



(&!)*(*) 



ji+rii 



ji! (ji + n iY- 



(78) 



where a*. n . = (— 1) 1 7^ \ The Hamiltonian in the 

interaction picture is given by V(t) = e lHot ^ h Hi nt e~ lHot ^ h , 
thus we have 



N 

(75) V(t) = M7cj + e- i(Adt+<w J| e - *"? H t (t) + H.c, (79) 



with 



+00 +00 



m)= E E 



rii = — 00 j i =ma,x[Q.-n i ] 



Ji! (ji + n,)! 



-iriiUJ-it 



(80) 
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where 



FIG. 7: f^,*. in Eq. d89t is plotted as the function of m for different 
values of r\ = 0.1, 0.3, 0.5, 0.7, 0.9 in (a) m = 0; (b) m = 1; (c) 
m — 10; and (d) m — 30. Here, for the convenience, we assume 
m = n and m, = m 



here A<j = cj^ — cjo- We assume that the off -resonant tran- 
sitions can be neglected in the resonant or near-resonant driv- 
ing condition. When the driving field satisfies the condition 

LO d = OJ - J^iLl We haVe 

N 

V {n * } = ftfio+e - ^" Y[ e-^H hrH + H.c, (81) 
with Hi >ni given by 



i=l 



ji+rii 



E 

ji— max[0, — ni] 



jil (ji +n l )l 



(82) 



With the Hamiltonian in Eq. dSTT i and using the complete- 
ness relation in Eq. d32l . the time evolution operator can be 
given as 



+oo 
{mj}=0 



(83) 



where 



Hereafter {n,} denotes a number series, that is, {rii} = 

ni, ri2, • • • , njv. The operator uf^X (t) can be further writ- 
ten as 



1 



c 



{mi- rii } 



1/2 



2\ 1/2 



1 



\e,{m,i}) (e,{rrii}\ 



(85) 



C 



with 



and the Rabi frequency 

njS } =nlla 

i=l 

If rii > 0, is given as 



(86) 



(87) 



(88) 



-if)? 

e 2 ? 



(m» + rii)! 



(^) 



(89) 



and if —rrii < rii < 0, we have 



S III , 



f_nl n *lf' n *l 



(90) 



here, L^. (77?) is the associated Laguerre polynomials. Be- 
cause is related to the effective Rabi frequency, thus it is 
plotted as a function of rii m Fig-LZl We find that decreases 
rapidly to the zero with the increase of m when rji <C 1 thus 
the Lamb-Dicke approximation is valid in this regime. How- 
ever, with the increase of rji, oscillates with the increase 
of rii and finally approach to the zero. Thus if rji is beyond the 
Lamb-Dicke regime, the terms for rii > 1 should be taken into 
account and these terms will make the state preparation much 
easier. It should be noted that £™\ decreases to zero within a 
finite number m even rji — 0.9. That is, there is maximum 
phonon number that we can create in one step even for the 
case that m is beyond the Lamb-Dicke regime. 

We now study the preparation of the NOON state of two 
mechanical resonators beyond the Lamb-Dicke regime. The 
time evolution operator of the system for two membranes in- 
side the cavity can be written out from Eq. (TS3T > by taking the 
subscripts and the superscripts {rrii} = mi, m^, {rrii ±rij} = 
mi ± iii,m2 ± ri2, and {rii} = ^1,^2- The coefficients in 
Eq. (l86l l for the case that two membranes inside the cavity are 



mi, 7712 



,-i<t>d 



(-1) 



ni+n 2 



\} L m u m 2 l ) 



with 



^mi ,7712 



/Vmi,n 2 
V mi ,m 



and the effective Rabi frequency 



077-1,77,2 _ 
7771 ,7772 



(91) 



(92) 



(93) 



Here, (™. is given by Eq. (|89l or Eq. (|90] l. We assume that 
the system is initially in the state \ip (to)) = \g, 0, 0). 
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In the step (i), the cavity field is driven by the external field 
which satisfies the carrier process. For the phase (j>d = 3ir/2 
and the time duration Aii = ir/2fl Q ), Q ) with the evolution op- 
erator U ' (Aii), the state of the system becomes 



\^(h)) = |e,0,0>. 



(94) 



In the step (ii), the frequency of the driving field is tuned to the 
TVth red sideband excitation corresponding to the frequency of 
the membrane one with = —Nlj%, With the time duration 
At2 = tt/AHqq and the phase <pd = 7r/2 + iV7r, the state of 
the system evolves to 



with the time evolution operator Uq'q'q (Aii). In the step (ii), 
the frequency of the driving field is tuned to the red sideband 
excitation such that A^ = — w-y — u)2 — ^3- With the time 
duration Ai2 = ^/^q'q'o an d the phase 4> c [ = 3tt/2, the 
state of the system evolves to 



|V»(t 2 )) = |fl)«»-^(|0,0,0) + |l,l,l». 



(101) 



with the time evolution operator U ' ' {At 2). Thus the sys- 
tem is prepared to a product state of the ground state \g) of the 
cavity field and the GHZ state of the three membranes. 



|^(t 2 )} = -^|e,O,O) + -^|0,i\r,O) 



(95) 



In the step (iii), the frequency of the driving field is tuned 
to the Nth red sideband excitation corresponding to the fre- 
quency of the membrane two such that A^ = —Nuj2- With 
the time duration A<3 = 7r/2f2Q and the phase <pd — 
tt/2 + Ntt, the state of the system evolves to 



1^3)} 



1 



\g)®—(\N,0) 



|0,iV». 



(96) 



Thus the NOON state of two mechanical modes is prepared 
with the cavity field in its ground state \g). It is obvious that 
the preparation process is more efficient than the one shown 
in the Lamb-Dicke regime. 

We now show the steps of generating the GHZ state. The 
time evolution operator of the three membranes inside the cav- 
ity can be easily given from Eq. ( f85T > by taking the subscripts 
and the superscripts as {m,} = mi, 1712,1713, {to,; ± m} = 
mi ± n\, ru2 ± 712, 7173 ± n^, and {rii} = ri\, 712, 713. The 
coefficients in Eq. (l86l l for the case of the three membranes 
are given as 



/^(ni,Tl2,l3 

mi,m2 ,7713 



-i<t>. 



fic-ir 



sin 



(r)n 1 ,n 2 ,n 3 
1 mi ) m2 ,1713 



(97) 



and 



with the Rabi frequency 



Qn 1 ,n 2 ,n 3 
mi ,m2,m3 



/iOrni,7i2,«3 Y 



(98) 



(99) 



Here, is given by Eq. d89l > or Eq. ( |90l >. 

We now assume that the system is initially prepared to the 
ground state \tp {to)) = \g, 0, 0, 0). In the step (i), the cavity 
field is driven to the carrier process. With the interaction time 
At 1 = ir / 4f2 ' ' and the phase <fid — 37r/2, the state of the 
system evolves to 



V2 



(|<?,0,0,0) + |e, 0,0,0)), 



(100) 



V. DISCUSSIONS AND CONCLUSIONS 

We discuss the experimental feasibility of the method by 
qualitatively considering environmental effects and informa- 
tion leakage, (i) The generation of entangled phonon states is 
based on the sideband excitations, which are extensively used 
in the optomechanical systems. Thus, our proposal should 
work in the resolved-sideband regime which requires every 
frequency Ui of the vibrational mode of mechanical mem- 
brane is bigger than the decay rate 7 C of the cavity field, 
i.e., Ui > j c . (ii) The two-level approximation in our pro- 
posal is guaranteed by the photon blockade effect, thus our 
method is more efficient when the single-photon strong cou- 
pling strength gi is much bigger than the decay rates j c and 
7m,i (i = 1) 2, • • • , N) of the cavity field and the mechani- 



cal modes, i.e., gi ^> j c , j m ,i- Also Eq. (1231 1 shows that the 
more number of the mechanical resonators corresponds to the 
higher nonlinearity of the cavity field, thus corresponds to the 
better two-level approximation of the cavity field, (iii) Dur- 
ing state preparation processes, negligible information leak- 
age from the ground or the first excited state to other up- 
per states of the cavity field requires that the anharmonic- 
ity ^2i{2gf /ijJi) of the cavity field induced by the mechan- 
ical modes should be much bigger than the strength f2 of 
the classical driving field in the carrier process 11751 . i.e., 
2^(2<?f /u)i) f2. (iv) To prevent information leakage due 
to nearly resonant transitions induced by different mechani- 
cal resonators, all of the transitions from the ground to the 
first excited state of the cavity field induced by the driving 
field and different mechanical resonators should be well sepa- 
rated in the frequency domain. In the the Lamb-Dicke approx- 
imation, the frequency differences between any two of mem- 
branes should satisfy the condition \u>i — u)j\ ^S> f2. However, 
beyond the Lamb-Dicke approximation as shown in Ref. H2I1 . 
the frequency differences between the membranes should sat- 
isfy the conditions |u;, 3> mm{u>i , ojj ) ^> fl ^> j c , 7^ 
or wia(wi,u}j) > |«j — Uj\ » > ToTm,*- From (iv), we 
can find that the membrane number that we can efficiently op- 
erate beyond the Lamb-Dicke approximation is much smaller 
than that in the Lamb-Dicke approximation. 

In summary, we have proposed a method to generate en- 
tangled states of vibrational modes of multiple membranes in- 
side the cavity via the radiation pressure. In particular, we 
carefully study the steps for generating several typical entan- 
gled states, e.g., the Bell and NOON state of two mechanical 
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(I) N is an odd number, (N-l)/2 is an odd number B 1 (R 1 B 1 ) [(N ~ 1)/2 ~ 11 ' 2 R 2 RiR 3 Ri— R N .iRiR N 

►|e,(N-l)/2,0,-,0> >- |g,0,0,-,0> ^ |g,0,-,0> 



|g,0,0, • • • ,0 > — >► |g,(N- 1 )/2,0, • • • ,0 >- 



C 



►|g,(N-l)/2,0,-,0> 



■> |e,N-l,0,— ,0> 



B^RiB,)^-'" 2 - 1 )' 2 



R2R1R3R1 "'Rn.iRiRn 



>► [&i,-,i> 
—\ — 



At, 



At, 



t 2 



At, 



At, 



t 4 Time 



(II) N is an even number, N/2 is an even number B^RiB,) 111 " 2 - 2 " 2 R 2 R 1 R 3 R 1 ---R N . 1 R 1 R N 

•|e,N/2-l,0,— ,0> ^ |g,0,0,-,0> ^. | g ,o,— ,0> 



|g,0,0, • • • ,0 > |e,N/2-l ,0, — ,0> ■ 



R, 



■|g,N/2,0,-,0> 



-> |e,N-l,0,-,0> 



Bi(RiB0 (N/2 - 2)/2 



R 2 RiR3Ri""R N .iRiR N 



> [g,l,"% 1> 



At, 



ti 



At, 



At, 



At 4 



Time 



(III) N is an odd number, (N-l)/2 is an even number (B 1 R,) <N - 1 > /4 R 2 R 1 R 3 R,-"R N . 1 R 1 R N 

►|e,(N-l)/2,0,-,0> >- |e,N-l,0,-,0> ^. |g,l,-,l> 

| g ,o,o,-,o>— >. | g ,(N-i)/2,o,-,o>-' 



•|g,(N-l)/2,0,-,0> 



> |g,0,0,-,0> 



> |g,o,-,o> 



(B,R,)" 



R2R1R3R] ""R]s[.iRiR>j 



At, 



At, 



At, 



At 4 



u 



Time 



(IV) N is an even number, N/2 is an odd number (BiR,) 1 * 2 -"' 2 R 2 RiR 3 R, •••R N _iRiR N 

' -|e,N/2,0 5 — ,0> ■ >- |e,N-l,0,—,0> >- [g,l,— ,1> 



|g,0,0,-,0> |g,N/2- 1,0,— ,0> - 



B, 



I g,N/2- 1,0, •••,()> 



> |g,0,0,-,0> 



(B,R,)' r 



R 2 RiR3Ri"-R N .iR 1 R N 



> |g,0,-,0> 



At, 



At, 



At, 



At 4 



Time 



FIG. 8: Schematic diagrams for preparing GHZ states of N membranes: (I) both N and (iV — l)/2 are odd numbers; (II) both iV and N/2 
are even numbers; (III) N is an odd number and (N — l)/2 is an even number; (IV) N is an even number and N/2 is an odd number. 



modes, the GHZ state and W state of three mechanical modes 
for the parameters with and without the Lamb-Dicke approx- 
imation. We should emphasize as following, (i) Basically, 
our method can be applied to other optomechanical systems in 
which many mechanical modes are coupled to a single-mode 
cavity field, such as optical cavity with levitating dielectric 
microspheres [96-99]] or tra pped atomic ensembles ll94l l95l . 
optomechanical crystals 1 30t 1 1 00^1 - and microwave cavity with 
nano-mechanical resonators M 8411 . (ii) Our proposal can be in 
principle used to produce any kind of the entangled states, 
(iii) We only qualitatively discuss the environmental effect 
and other information leakage on the state preparation. The 



quantitative analysis on these factors will be given in else- 
where, (iv) Our proposal is experimentally possible when 
the optomechanical coupling strength approaches the single- 
photon strong coupling limit. 
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Appendix A: Preparation of W and GHZ states of TV can be generated by sequentially applying a series of red- 

membranes sideband excitations, blue-sideband excitations and carrier 

Follow the method we have given in Sec. Ill, the W and 
GHZ states of N membranes inside the cavity as 

1 „, 
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line operations. In this appendix, we are going to give the 
method for generating the W and GHZ states of N membranes 
inside the cavity under the Lamb-Dicke approximation. 

Assume that the system for N membranes in the cavity is 
prepared in the state \if) (to)) = \g, 0® N ) initially, then the 
W state of the N membranes can be generated by one carrier 
process followed by N red-sideband excitations. After the 
action of the carrier process for time duration Ati = it/2Q,, 
the state of the system becomes 



\1>(t 1 )) = \e,0* ri ) 



(A3) 



at time t\ = to + Ati. After that, N red-sideband excita- 
tions (driving fields) are applied to the system sequentially, 
with frequencies ujj = lu — Wj, time durations Ai i+1 = 
[arcsin (l/^/N + 1 - 

1, • • ■ , N. At time t 2 
system becomes 



an d phase 



7r/2 for i 



t\ + Yli=i A^i+i, the state of the 



1 



\^(t 2 )) = | 5 )®-=(|l,0,0,---,0) + |0,l,0,---,0) 
V A 



+ ••• + 10, 0,--- ,0,1)). 



(A4) 



As shown in Fig. [8] we can also generate the GHZ state of 
N membranes for the system is initially prepared in the state 
\4> (to)) = \g,0® N ). In step (i), by sequentially applying a 
series of red-sideband and carrier processes, we can prepare 
the state 



\i>(h)} 



N - 1 



0®(JV- 



for the odd number N, or 



liK*i)> 



e, 



N 



1 ,o l 



for the even numbers N and N/2, or 



1-0 (*i)> = 



//•(f -i).o iX 1 



(A5) 



(A6) 



(A7) 



for the even number N and odd number N/2. 

In step (ii), if N is an odd number, the system can be pre- 
pared in the state 



1 



|^ 2 )) = _ (Ig) + |e)) ® 



N - 1 



(A8) 



by the action of the carrier process for time duration At 2 = 
tt/AH; if both N and N/2 are even numbers, the system can 
be prepared in the state 



hH*2)> = 



i 

V2 



N 



N 



9, 



i(N-l) 



(A9) 

by the red sideband excitation of the membrane one for time 
duration At 2 = tt / 4£l^ v ^ 2 , and if N is an even number but 
N/2 is an odd number, the system can be prepared in the state 



hH*2)> 



V2 



N 



N 



(A10) 




10 N 15 



FIG. 9: Fidelity for preparing GHZ state -±= (|0® w ) + ll®^)) 
from N = 3 to N = 22. 



by the blue sideband excitation of the membrane one for time 
duration At 2 = 7r/4f2^y 2 . 

In step (iii), we can prepare the state 



|V> (t 3 )) =-j=(\ g) Q) + \e,N-l)) 



(All) 



for all the cases in step (ii) by a series of processes as shown 
in Fig. [8] There are some information leaked in this step for 
the states \g, 0, 0®^" 1 )) and \e,N- 1, 0®^ -1 )) can not be 
prepared in synchronism, and the fidelity is shown in Fig. [9] 
(discussion see the Appendix B). 

In step (iv), by the action of R 2 R\R^Ri ■ ■ ■ Rn-\R\Rn 

for time duration At± = X^? 1 7r /2^i + Yli= 2 7r/20|' we 
get the state 



\^(t 4 )} = \g)® — (\0t 



r >+|l® W » 



(A12) 



Thus we have prepared N membranes in the GHZ state and 
the optical field in the cavity in the ground state \g). 



Appendix B: The information leakage caused by 
non-synchronization 

In the preparation of the NOON state and GHZ state in 
the Lamb-Dicke regime, we need the processes for prepa- 
ration two different states synchronic ally. For example, in 
the preparation of state (|2,0) + 1 , 2) ) /v2, the transitions 
\g, 1, 1) -> |e, 1, 0), and |e, 0, 1) -> \g, 0, 2) should be syn- 
chronized. In the other words, the time duration should sat- 
isfy Eq. (l54l . However, we find that Eq. (l54l i can only be 
satisfied in the infinite approximation. In this appendix, we 
will discuss the information leakage caused by this type of 
non-synchronization. 

Without lose of generality, suppose that the state of the sys- 
tem at time t is given as 

|^ (t)) = ±= |e, n - 1, n') + -±= \g, m, m') (Bl) 
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a (n',m') 




10 , 15 20 



FIG. 10: The maximal value of a versus n and m . 



and at the following time t' , we need to prepare the state 

|V (i')> = |g, n, «0 + |e, m - 1, m') . (B2) 

So we drive the optical cavity with a laser field resonant to 
the first red sideband of the first membrane for time duration 
At = t' — t, and the system evolves into 

hK*0) = -^[cos(fiiAt)|e,n-l,n') 

-ijti sin (ft* At) |p,7i,n') 

+ cos (fi^Ai) \g,m,m') 

-ie-^r sin (fi^At) |e, m - 1, to')]- (B3) 

where ft* = ft?7i -y/n, ft^ = Qr]i \fm. 

In order to make sure that the system at the time t' is in the 
state given by Eq. dB2l) . the time duration At should satisfy the 
equations 

sin (ft* At) = ±1, sin (ft*„At) = ±1. (B4) 
Eq.( IB4l ) can be equal to the equations 

(n 1 n ±nl n )At = p7r,n 1 m At = (q + ±)ir, (B5) 

where p and q are positive integral numbers. Taking 



nl m 

ten as 



yjn/m (suppose n > to), Eq.( IB5l ) can be rewrit- 



± 1 



P- 



(B6) 



for choosing appreciate integral numbers p and q. 

According to the rate of the two angular frequencies, the 
solutions can be divided into three cases: 



(i) If y/n/m is a rational number, we can rewrite it as a 
quotient of integers, 



■ — (B7) 

TO TO' 

where n' and to' are integers and they have no common fac- 
tors. Substitute Eq. (|B7| > into Eq. (|B6K we get 



n'±m'\ {2q+l 



(B8) 



If both n' and to' are odd, then there are integral numbers p 
and q that satisfy Eq.( |B8l ) exactly. 



(ii) If y/n/m = n'/m' is a rational number, and one of 
them (n' and m') is even, there are no integral numbers satis- 
fying Eq. dB8l ). However, there are integral numbers can sat- 
isfy the following equations for appreciate value of At, 



sin (ft* At) = ± sin (ft^ At) = ±a, 



where a is given by 



sin 



m 



n' ± to 



(B9) 



(BIO) 



By the increase of p, a oscillates periodically with a period 
n' ± to'. The maximal value of a as function of n' and m' (de- 
noted by a max (n', to')) is shown in FigJlO] From the figure 
we can see that a max (n', to') > 0.86, and a max (n' , to') = 1 
for the case that both n' and m' are odd, which is agree with 
the result given above. 

(iii) If y/n/m is an irrational number, there are no integral 
numbers (p and q) that satisfy Eq. (IB6l ) exactly, but there is a 
rational number n'/m! (both n' and to' are odd) that can be 
infinitely close to y/n/m. So we can get appreciate value of 
At which can satisfy Eq.( IB6l ) in the infinite approximation. 
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